Chaos and statistical relaxation in quantum systems of interacting particles 
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We propose a method to study the transition to chaos in isolated quantum systems of interacting 
particles. It is based on the concept of delocalization of eigenstates in the energy shell, controlled 
by the Gaussian form of the strength function. We show that although the fluctuations of energy 
levels in integrable and non-integrable systems are principally different, global properties of the 
eigenstates may be quite similar, provided the interaction between particles exceeds some critical 
value. In this case the quench dynamics can be described analytically, demonstrating the universal 
statistical relaxation of the systems irrespectively of whether they are integrable or not. 
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Introduction. - Recent experimental progresses in the 
studies of various quantum systems of interacting par- 
ticles (see, e.g., [H) have triggered the interest in basic 
problems of many-body physics. One of the issues that 
has been widely discussed in the literature is the onset of 
thermalization in an isolated quantum system caused by 
interparticle interactions |2h10I]. 

A prerequisite for thermalization is the statistical re- 
laxation of the system to some kind of equilibrium and 
its viability has been associated with the onset of quan- 
tum chaos. The latter term was originally created to ad- 
dress specific properties of dynamical quantum systems 
whose classical counterparts are chaotic. Later, it was 
found that similar properties of spectra, eigenstates and 
dynamics could emerge in quantum systems without a 
classical limit, as well as in quantum systems with disor- 
dered potentials. Nowadays, the term quantum chaos is 
used in a broader context when referring to those prop- 
erties, irrespectively of the existence of a classical limit. 

According to studies of isolated quantum many-body 
systems, their eigenfunctions (EFs) in the mean-field ba- 
sis spread as the interaction between particles (or quasi- 
particles) increases, and they may eventually become 
chaotic eigenstates [11,]. A crucial point is that due to 
the finite range of the interactions, only part of the un- 
perturbed basis states | n) is directly coupled and therefore 
accessible to the eigenstates |a} — J2n I"-)- In the en- 
ergy representation, this fraction constitutes the energy 
shell of the system, being either partly or fully filled by 
the actual eigenstates [12|. In the first case, the eigen- 
states are localized, having a small number of non-zero 
elements C^. In the second case, when the energy shell 
is full y fi lled, the eigenstates can be either sparse or er- 
godic [l3| , both showing a very large number of principal 
components {Npc >> 1) strongly fiuctuating with n. In 
ergodic eigenstates the coefficients C" become random 
variables following a Gaussian distribution around the 
"envelope" defined by the energy shell. The latter sce- 



nario is that of chaotic eigenstates and occurs when the 
interaction exceeds a critical value 12 1. 

As argued in Ref. 12|, the energy shell is associated 



with the limiting form of the strength function (SF) for 
a given interaction strength. This function is defined 
via the projection of unperturbed states onto the basis 
of perturbed (exact) eigenstates. Written in the energy 
representation, SF is widely used in nuclear physics and is 
analogous to local density of states in solid state physics. 
It has also been shown (see, e.g., 0, d, [l3, US] ) that the 
shape of SF changes from Breit-Wigner (Lorentzian) to 
Gaussian as the interparticle interaction increases. 

It should be stressed that if a quantum system has a 
classical limit, the shapes of both EFs and SFs in the en- 
ergy representation have classical analogs. The quantum- 
classical correspondence of EFs and SFs has been studied 
in great detail for various few and many-body systems 
(see Ref. Typically, delocalization of EFs in the 

energy shell is directly related to the chaotization of the 
system in the classical limit, thus providing a tool to re- 
veal the transition to quantum chaos. 

The emergence of chaotic eigenstates has been related 
to the onset of thermalization in many-body systems, 
even if the latter are isolated from a heat bath [2t9| • An 
essential point of Refs. 0-0] is that when the eigenstates 
become chaotic, the distribution of occupation numbers 
achieves a Fermi-Dirac or a Bose-Einstein form, thus al- 
lowing for the introduction of temperature. Using a two- 
body random matrix model, a relation between tempera- 
ture and interaction strength was analytically derived @ . 
Therefore, the interparticle interaction plays the role of 
a heat bath when the eigenstates are chaotic. Since the 
components of chaotic eigenstates can be treated as ran- 
dom variables, the eigenstates close in energy are statisti- 
cally similar. This fact has been employed to justify the 
agreement between the expectation values of few-body 
observables and the predictions from a microcanonical 
ensemble 
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In this Letter we present a unified approach to iden- 
tify the emergence of chaotic eigenstates. Two models of 
interacting spins-1/2 are considered: one model is com- 
pletely integrable for any value of the interaction and 
the other is not. The different procedures used to deter- 
mine the critical parameters above which the eigenstates 
become chaotic lead all to the same values. Our results 
reveal a relation between energy level statistics, the struc- 
ture of EFs and SFs, and the quench dynamics. We show 
that statistical relaxation occurs for both models and can 
be described analytically with the same expression. 

The models. We consider two models of interact- 
ing spins-1/2. One model has only nearest- neighbor 
couplings, which results in the complete integrability 
of the system. The other case has additional next- 
nearest-neighbor couplings, and becomes chaotic when 
the two coupling strengths are comparable. Assuming 
open boundary conditions, the Hamiltonians read as 

Hi^Ho + nVu H2 = Hi+XV2, (1) 

L-l L-1 

Ho = {SfSf_^i + SfSf^j^) , Vi ^ S-S^^i, 
1=1 1=1 

L-2 
i=l 

where fj, and A control the perturbation in Model 1 and 
Model 2, respectively. Here, L is the number of sites, 
g^'V'^ — cr^'V'^ /2 are the spin operators at site i, with 
tjf'*''^ as the Pauli matrices and ?i = 1. The coupling 
strength J defines the energy scale and is set to 1. We 
refer to a spin "up" in the z direction as an excitation, 
assuming the presence of a magnetic field. 

In Model 1, the term Hq determines the (mean field) 
basis in which the total Hamiltonian Hi is presented. 
This term moves the excitations through the chain and 
can be mapped onto a system of noninteracting spinless 
fermions [16[ or of hard-core bosons 17 1 , being therefore 
integrable. The system remains integrable even when 
the Ising interaction Vi is added, no matter how large 
the anisotropy parameter ^ is. The total Hamiltonian 
Hi is known as the XXZ Hamiltonian and can be solved 



exactly via the Bethe Ansatz [18 1. 

The integrable XXZ Hamiltonian determines the mean 
field basis for Model 2. Thus, V2 is treated as the "resid- 
ual interaction" responsible for the onset of chaos. The 
parameter A refers to the ratio between the next-nearest- 
neighbor and the nearest-neighbor exchange. 

Depending on the parameters of the Hamiltonians ([1} , 
different symmetries are identified [19,]. For the sake of 
generality, we avoid them by restricting our analysis to 
a subspace with L/3 up-spins and /i 7^ 1. The only re- 
maining symmetry is parity. We take it into account by 
studying only even states, which leads to subspaces of 
dimension TV - (1/2)L!/[(L/3)!(L - L/3)!]. All data are 
given for L — 15 and /i = 0.5 for Model 2. 



Spectrum statistics. - According to the common lore, 
we analyze first the level spacing distribution P{s) for 
both models, numerically obtained for different values of 
the control parameters /i and A. For Model 1, P{s) is 
close to the Poisson distribution for any value of /i. For 
Model 2, the transition of P{s) from Poisson to Wigner- 
Dyson as A increases is shown in Fig. [TJ The standard 
approach of fitting P(s) with the Brody distribution [2^ 
allows us to extract the repulsion parameter /3 character- 
izing the transition between the two distributions. From 
Fig. [I] the transition for Model 2 with L = 15 occurs 
approximately at A ~ 0.5. 




FIG. 1: Left: P{a) for Model 2 with A = 0.1, 0.5 compared 
with the Wigner-Dyson distribution (smooth curve). Right: 
Brody parameter /3 as a function of A. 

Emergence of chaotic eigenstates. - Much more infor- 
mation is contained in the structure of the eigenstates. 
Our data show that as the strength of the perturbations 
Vi and V2 increases, the eigenstates of both integrable 
and non-integrable models undergo a transition from lo- 
calized to chaotic-like ones. Typical examples of the am- 
plitudes C" of such eigenstates with energy Ea from the 
center of the energy band are shown in Fig. [2l Here, the 
eigenstates are given as a function of the unperturbed 
energy e„ rather than in the basis representation, follow- 
ing the one-to-one correspondence between each unper- 
turbed state \n) and its energy The figure shows how 
the eigenstates spread as the interaction increases. 
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FIG. 2: Typical localized (top) and extended (bottom) eigen- 
states for Model 1 (left) and Model 2 (right). 

In order to find the critical parameters /icr and Xcr 
above which the perturbation is strong and the eigen- 
states are extended (in energy shell), different approaches 
may be employed. We start by analyzing the matrix ele- 
ments of Hi and H2. It is important to take into account 
that in each line n of the Hamiltonians, the perturbation 
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couples directly only M„ unperturbed states 0, [2lj . We 
have numerically found that at the center of the energy 
band, M„ w N/2 and ^ N for Models 1 and 2, 
respectively. Thus, the Hamiltonian matrix of the inte- 
grable model has more vanishing elements i?„m, being 
more sparse than the matrix of Model 2. 

To determine the critical perturbation, we compare 
the average value of the coupling strength, Vn — 
^m^n \Hnm\lMn, of each Hue n with the mean level 
spacing c?„ between directly coupled states. The mean 
level spacing can be estimated as d„ = [e™""^ — e™*"]/M„, 
where e™"^ (e™™) is the unperturbed energy correspond- 
ing to the largest (smallest) m where Hnm 7^ 0. Our 
results show that for ^ > ^cr ~ 0.5 and A > Xcr ~ 0.5 
the ratio Vn/dn becomes larger than 1 and the pertur- 
bation is considered to be strong. Notice that for Model 
2, the obtained value of Xcr corresponds to the onset 
of the Wigner-Dyson statistics, as independently found 
from the level spacing distribution. This is remarkable if 
we take into account that no diagonalization was neces- 
sary to derive the above estimates. 

Strength function: From Breit- Wigner to Gaussian. - 
Another way to obtain the critical values /icr and Xcr re- 
lies on the shape of SF. The latter corresponds to the 
dependence of — |C"p on the exact energies Ea 
for each fixed unperturbed energy £„. It contains in- 
formation about the energies Ea that become accessible 
to an initial state \n) when the perturbation is turned 
on. Clearly, SF is related also to the structure of EFs, 
since the latter is derived from the same , but now as 
a function of the unperturbed energies e„. 

In quantum many-body systems, the form of SF typi- 
cally changes from Breit- Wigner to Gaussian as the inter- 
particle interactions increase 0, [^, [13] ■ This transition 
occurs when the half-width of the Breit- Wigner distri- 
bution becomes comparable to the width of the energy 
shell. In this case, as we show next, there emerge chaotic 
eigenstates filling the whole available energy shell. 

The energy shell corresponds to the distribution of 
states obtained from a matrix filled only with the off- 
diagonal elements of the perturbation. It is associ- 
ated with the maximal SF, that is the SF that arises 
when the diagonal part of the Hamiltonian can be ne- 
glected. We verified that the energy shell coincides with 
a Gaussian of variance ct^ given by the second moment 
of the off-diagonal elements of the matrix Hamiltonian, 
'^^ — 'Ylira+n \^nrn^ Note that uo diagonahzation is 
required to derive this expression. 

Our numerical data confirm that the transition from 
Breit- Wigner to Gaussian occurs for the same critical val- 
ues obtained above, \icr , Xcr ~ 0.5, as indicated in Fig.[3l 
In the figure, the envelopes of the SFs were obtained by 
smoothing the dependence of w" on Ea for fixed unper- 
turbed energies £„ with n « N/2. The fit to either a 
Breit- Wigner or a Gaussian form was done with high ac- 
curacy, allowing us to discriminate between the two func- 




FIG. 3: (Color online) Strength functions for Model 1 (left) 
and Model 2 (right) obtained by averaging over 5 states in 
the middle of the spectrum. Circles give a Breit- Wigner fit 
(fi = X — 0.4) and a Gaussian fit {fj, — 1.5 and A — 1.0). Solid 
curves correspond to the Gaussian form of the energy shells. 



tions. It is noteworthy the excellent agreement between 
the Gaussian fit and the Gaussian obtained simply from 
the off-diagonal elements of the Hamiltonians. 

Structure of eigenstates in energy shell. - The eigen- 
states may be localized, sparse or ergodically extended 
in the energy shell. The data in Fig. 0] demonstrate that 
for a sufficiently strong perturbation the eigenstates un- 
dergo a transition from strongly localized to extended 
states, somehow filling the energy shell. The transition 
to chaotic eigenstates occurs again at the same critical 
parameters iicr , Xcr ~ 0.5. These results confirm the 
predictions made on the basis of both, the estimate of 
Vn/dn and the Gaussian form of the strength functions. 
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FIG. 4: Structure of eigenstates in the energy shells for Model 
1 (left) and Model 2 (right) obtained by averaging over 5 
states in the middle of the energy band. Solid curves corre- 
spond to the Gaussian form of the energy shells. 

One notices that above the critical value, the eigen- 
states of Model 2 appear to fill the energy shell better 
than those of Model 1 (cf. bottom panels of Fig. U]). 
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Indeed, for Model 1, the eigenstates are not completely 
extended in the energy shell even for strong interaction, 
^ = 1.5, although they do fill a large part of it. At the 
same token, a close inspection of the level of delocaliza- 
tion of individual EFs and SFs has revealed differences 
between the two models. Overall, delocalization mea- 
sures, such as the inverse participation ratio or Shannon 
entropy, show larger fluctuations for Model 1 than for 
Model 2 [22]. This agrees with recent results obtained 
for bosonic and fermionic systems 

Statistical relaxation. - The knowledge of the shape 
of SFs allows one to describe the quench dynamics in a 
system of interacting particles, provided the eigenstates 
are chaotic and delocalized in the energy shell j23|]. By 
quench dynamics we mean the time evolution of initial 
states corresponding to unperturbed basis vectors which 
takes place once the interaction is turned on. To see how 
the relaxation occurs, we study the time dependence of 
the Shannon entropy S in the unperturbed basis. For an 
initial state \k), it is defined as 



N 



s,{t)^-J2^i!:Ht)^nni'Ht) 



(2) 
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FIG. 5: (Color online) Shannon entropy vs rescaled time for 
Model 1 (left) and Model 2 (right) for strong perturbation. 
Circles stand for numerical data, solid curves correspond to 
Eq. (131), and dashed lines show the linear dependence 

An analytical expression for S{t) has been derived with 
the use of a two-body random matrix model ^2^ , 



Skit) = ~Wk{t) \nWk{t) - [1 - Wk{t)] In 



Wk{t) \ 



pc 



(3) 

Here Wk (t) is the probability for the system to stay in the 
initial state \k) and Npc is the average number of directly 
coupled states. We obtain Npc numerically according to 
^pc = (e'^), where the average (.) is performed over a 
long time after the entropy saturates. As for Wk{t), it 
has been shown [ij] that for a Gaussian SF, it decays as 
Wk{t) = exp(— cr^t^), with cr^ as previously defined. 

Figure [5] shows numerical data for the relaxation pro- 
cess of both models. To reduce fluctuations, we aver- 
age over 5 initial basis states excited in a narrow energy 
range in the middle of the spectrum. Initially, the en- 
tropy grows quadratically, as given by perturbation the- 
ory. Afterwards, a clear linear growth is observed before 



S{t) reaches relaxation. With high accuracy the linear 
behavior of S{t) is described by the simple relation [23{, 



Sk{t)^aktlnMk. 



(4) 



Note that Eq. Q depends only on the elements of the 



Hamiltonian: az 



and Mk is the con- 



nectivity, that is the number of directly coupled unper- 
turbed states in the fc-th line of the Hamiltonian matrix. 
As seen in Fig. [SJ the analytical expressions ^ and @ 
give a correct description of the increase of the entropy 
for both models in the regime corresponding to the onset 
of chaotic-like eigenstates delocalized in the energy shell. 
The same relation Q was found to emerge also for an 
integrable model of interacting bosons 2J| . 

Conclusion. - We have studied spectrum statistics, 
structures of eigenstates and strength functions for two 
models of interacting spins, connecting the results with 
the onset of chaotic-like eigenstates and the statistical 
relaxation of the quench dynamics. The key point of 
our approach is the existence of an energy shell of finite 
range, inside which the eigenfunctions can be either lo- 
calized or extended. We have shown that the critical pa- 
rameters above which the eigenstates become chaotic can 
be equally found by simply studying the elements of the 
Hamiltonian or by analyzing the shape of the strength 
function. The latter provides us with the form of the en- 
ergy shell, thus allowing one to clearly define the notion 
of delocalized eigenstates in the energy shell. 

By studying the time dependence of the Shannon en- 
tropy, we have shown that numerical data are in full 
agreement with the analytical predictions of the quench 
dynamics, provided the eigenstates are chaotic-like. Our 
approach is very general and expected to apply to dif- 
ferent systems of interacting particles, such as those cur- 
rently under theoretical and experimental investigation. 
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